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Abstract This study based on longitudinal Zeeman effect magnetograms and 
spectral line scans investigates the dependence of solar surface magnetic fields 
on the spectral line used and the way the line is sampled in order to estimate 
the magnetic flux emerging above the solar atmosphere and penetrating to the 
corona from magnetograms of the Mt. Wilson 150-foot tower synoptic program 
(MWO). We have compared the synoptic program A5250A line of Fel to the 
line of Fe I at A5233A since this latter line has a broad shape with a profile that 
is nearly linear over a large portion of its wings. The present study uses five 
pairs of sampling points on the A5233A line. Line profile observations show that 
the determination of the field strength from the Stokes V parameter or from 
line bisectors in the circularly polarized line profiles lead to similar dependencies 
on the spectral sampling of the lines with the bisector method being the less 
sensitive. We recommend adoption of the field determined with the line bisector 
method as the best estimate of the emergent photospheric flux and further 
recommend the use of a sampling point as close to the line core as is practical. 
The combination of the line profile measurements and the cross-correlation of 
fields measured simultaneously with A5250A and A5233A yields a formula for the 
scale factor 5^^ that multiplies the MWO synoptic magnetic fields. Using p as the 
centcr-to-limb angle (CLA), a fit to this scale factor is ~ 4.15 — 2.82 sm^{p). 
Previously 6'-^ = 4.5 — 2.5 sin^(p) had been used. The new calibration shows that 
magnetic fields measured by the MDI system on the SOHO spacecraft are equal 
to 0.619 ± 0.018 times the true value at a center-to-limb position 30°. Berger 
and Lites (2003) found this factor to be 0.64 ± 0.013 based on a comparison the 
the Advanced Stokes Polarimeter. 

Keywords: Chromosphere, Active; Magnetic fields, Chromosphere; Magnetic 
fields. Photosphere 

1. Introduction 

The sun's 22-year cycle of activity is most clearly seen in solar surface magnetic 
fields (see for example the review by Ossendrijver (2003) for a good discussion 
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of the solar dynamo and solar magnetic fields). Although the Sunspot Number 
(SSN) is the most readily available indicator of the state of the solar cycle, the 
magnetic fields over the whole surface provide a more complete measure of the 
dynamo process. A variety of questions arise in using the surface magnetic fields 
to study the solar cycle - which fields are most important, the weak general field 
or the strong field and associated sunspots, how do the field strengths change 
over various time scales from instabilities with changes in seconds or less to 
trends lasting many decades. The connection between photospheric magnetic 
fields and magnetic fields near earth or at interplanetary spacecraft requires 
knowledge of the field strength at the solar surface. Is the overall strength of 
the magnetic field stationary or does it have any multi-decade trends (cf: Arge 
et al. (2002))? These surface fields are also needed to calculate the strength 
measured by interplanetary spacecraft such as Ulysses (Gibson et al., 1999). The 
topology of the field in the solar surface regions introduces additional uncertainty 
since the field of one polarity lower in the atmosphere may be canceled by the 
opposite polarity before it can emerge into the corona and hcliospherc (Antiochos 
et al., 2007). Questions such as these involving magnetic field strength can only 
be addressed when we are confident in the quantitative interpretation of the 
measured quantity. 

The f 50- foot solar tower telescope on Mt. Wilson has been dedicated to the 
synoptic measurement of solar surface magnetic fields using a Babcock magneto- 
graph observing the Fe I spectral line at A5250A. Although other spectral lines 
and observing techniques offer advantages over these choices, the uniformity of 
the data set extending back to 1967 is a strong reason to leave the observing 
system as nearly unaltered as possible. The application of these data to modeling 
the heliospheric magnetic field then depends on our ability to understand how 
to interpret the magnetograms. An early step toward this goal was carried out 
by Howard and Stcnflo (f 972) and subsequently extended by Ulrich (1992) and 
Ulrich et al. (2002) using a cross-correlation of fields observed simultaneously 
from two or more pairs of spectral samples. A compelling feature of the scatter 
diagrams from these simultaneous measurements is the linearity between the 
magnetic fields obtained from the different spectral samplings. Each spectral 
sample measures a different altitude in the solar atmosphere and is influenced 
by different thermal properties of the atmosphere at that altitude. If the altitude 
dependent part of the relationship were also to depend on the field strength at 
each point on the solar surface, the scatter diagram would not be confined to 
a straight-line band passing through the zero-zero point. The observed straight 
line distributions found on the scatter diagram establish that variable filling 
factor is the only cause for varying field strength, all flux tubes being essentially 
identical. The fact that the slope is not unity can be used to study the altitude 
dependence of the flux tubes. 

Throughout this paper we will use the term "magnetic held strength" and 
usually mean the apparent field strength which is the product of a flux tube 
field strength and a filling factor. This quantity is actually the magnetic fiux 
per unit area and is the quantity we wish to measure. The extensive literature 
devoted to the determination of the flux tube fields and the filling factors is of 
interest to us and is referenced where appropriate below. However, questions of 
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flux tube physics are of secondary importance to our objectives in this paper. 
Where the flux tube field strength is discussed, it is made explicit. Otherwise, 
the term magnetic field strength should be read to mean the magnetic flux per 
unit area as measured in the observed pixels. 

The essential step in our study of magnetic fleld strength is the choice of a 
comparison line. Our objective is to learn the relationship between measured 
circular polarization of a relatively large pixel (12 arcseconds squared or 20 
arcseconds squared) and the magnetic flux in that area. We are less concerned 
with questions of the structure of magnetic flux tubes or the strength of the field 
within the flux tubes - we only seek to learn the best way to estimate the product 
of the flUing factor and the magnetic field in the fiux tubes. In addition, we need 
to be able to observe the line simultaneously with A5250A. For this purpose the 
line at A5233A is ideally suited due to its large width and depth and due to 
the good linearity of its line wings. For a related objective of determining the 
flux tube magnetic fleld strength the comparison line chosen by Stenflo (1973) 
was A5247A to take advantage of the fact that the two lines have differing g^s 
values but very similar thermal response and line formation properties. However 
in both this paper and the subsequent study by Frazier and Stenflo (1978) the 
reference to A5233A was retained with the assumption that magnetic fluxes from 
this line are obtained directly from the polarization factors and the line proflle. 
An interesting citation to a private communication by Livingston in this last 
paper suggests that the magnetic fields derived from A5233A depend on position 
within this line. We verify below that the deduced magnetic field strength for 
any particular area on the sun's surface that depends on the spectral sampling 
of A5233. Our previous results based on A5233A yielded a correction factor of 
4.5 to be applied as a multiplier for the magnetic fields obtained from A5250A 
whereas Stenfio (1973) found only a factor of two. Our correction factor is larger 
due partly to the fact that the shifted Zeeman components at A5250A are shifted 
beyond the sampling pass-band of the MWO Babcock magnetograph. We refer to 
the reduced apparent field due to the large line shift as the "classical saturation 
effect" since it depends only on the actual line profile, the strength of the flux 
tube magnetic fields and the spectral resolution of the observing system. 

2. Magnetogram Observations 

As a first step in establishing a more complete understanding of the scale factors 
relating different observing methodologies we have carried out a new series of 
observations utilizing simultaneous magnetic field measurements with six differ- 
ent spectral configurations. The MWO system integrates all spectral channels 
continuously with all integrator registers being emptied sequentially at a rate of 
40 Hz. Due to this feature of the system image motion and smearing affects all 
spectral samples identically so that the spatial average of spatial inhomogenieties 
is also identical. Thus ratios of magnetic fields calculated from different sampling 
pairs yield relative magnetic fields from identical solar surface structures. Start- 
ing with the concept that the magnetic field from A5233A indicates the true 
magnetic fiux, the ratio ~ ^a5233A/^A525oA should be the scale factor that 
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converts the MWO synoptic program magnetic fields to the appropriate value 
emerging from the photospheric layers. For the actual observations we need to 
specify the spectral sampling as well as the spectral line. The study by Ulrich 
(1992) used A5250A±39mA and A5233A±45mA to yield the scale factor of 4.5. 
Contrary to the assumptions of Frazier and Stenflo (1972), Stenflo (1973) and 
Frazier and Stenflo (1978) we find that the field measured by A5233A depends 
on the spectral sampling so that we must generalize the notation. We do this 
by defining a new quantity f?J^'^ which relates the magnetic field for the spectral 
line at wavelength a having spectral sampling offset by n to the magnetic field 
for the spectral line at wavelength b having spectral sampling of m with a and 
b in Angstroms and n and m in niA. Thus we write: 

O fl-^i 71) 

so that the previously defined is: 

r-l _ 5233,45 
" ~ %250,39 • 

The fact that the magnetic field derived from A5233A depends on the sampling 
separation is not a surprise in view of the complexity of the atmospheric flux 
tube structure but it does present us with the additional task of determining 
which portion of this line is most appropriate as an indicator of the emergent 
magnetic flux. 

During the interval 18 April 2007 to 8 May 2007 a set of five special magne- 
tograms were obtained using the 24-channel system described by Ulrich et al. 
(2002) with the blue, 10-channcl fiberoptic reformattor centered on A5233A 
instead of the A5896A line of NaD. This configuration provided us with five line 
pairs having separations of ±9, ±29, ±84, ±102 and ± 177mA. The observations 
in 1992 were made using a different spectrograph configuration wherein the 
spectral line sampling was at ±45mA. The sampling of A5250A was left as in 
the synoptic program at ±39mA. During the scanning process the light entering 
the spectrograph is unfiltered. The dispersed light passes through narrow-band 
blocking filters where it falls on all spectral sampling pickups simultaneously. 
The output currents from all photomultiplier tubes are summed individually and 
simultaneously into a set of registers. We refer to the sequence from the entry 
window of a single spectral sample to the summing register as a channel. The final 
digitization numbers for all channels for each pixel then originate from identical 
portions of the solar surface even though the sampled area is not precisely defined 
by the entrance aperture due to distortions of the solar image by the Earth's 
atmosphere. We compute magnetic fields from each spectral sampling pair and 
compare the derived fields in a scatter-diagram format. Since we expect the 
interrelationships to depend on position within the solar image, we restrict the 
scatter diagrams to pixels having a specified range in center-to-limb angle p. 

The comparisons were done between the various samplings and the one at 
A5233A ± 84mA. The corresponding scatter diagrams for three of the five cases 
are shown in Figures 1 to 3. The values of ^75233'™ are given in Table 1 for all five 
cases. The adoption of n = 84 is somewhat arbitrary but was done because the 
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Figure 1. This figure shows the relationship between magnetic fields determined from spectral 

samples at A5233A ± 29mA and at A5233A it 84mA. Each point on the figure shown as the 

very small X represents a single pixel observed once. The range in ccntcr-to-limb angle p for 

each figure is shown above each figure. The light at the sampling wavelengths is recorded over 

identical time and space limits so that the difference in determined field strengths is a result 

of differing physical conditions at the two different heights of formation. The portion of the 

figure near the 0,0 origin is without points in order to improve the manageability of the plot by 

decreasing the number of points plotted. The slope of the line relating the two determined field 

5233 84 

strengths is indicated in the upper left section of each plot. This quantity '(5233 '29 multiplies 
the field determined from A5233A ± 29mA in order to convert that field into the value that 
would be determined from the combination at A5233A ± 84mA. The scatter of the points is 
due to photon statistics as is discussed in the text. 
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Figure 2. This figure is in the same format as Figure 1 but shows the relationship between 
magnetie fields determined from spectral samples at A5233A± 102mA and at A5233A±84mA. 
The tightness of the relationship is noteworthy. 



Table 1. The correlation slopes V5233'm ^^'^ their errors as a function of m and sin(p). 



m 


< sin{p) < 0.4 


0.4 < sin(p) < 0.6 


0.6 < sin{p) < 0.8 


0.8 < sin{p) < 0.95 


9 


1.665±0.067 


1.887±0.043 


2.262±0.077 


2.042±0.105 


29 


1.795±0.027 


1.836±0.030 


2.169±0.041 


2.101±0.053 


102 


1.162±0.010 


1.133±0.012 


1.134±0.010 


1.009±0.019 


177 


2.694±0.094 


2.727±0.085 


2.706±0.075 


2.607±0.206 
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Figure 3. This figure is in the same format as Figure 1 but shows the relationship between 
magnetic fields determined from spectral samples at A5250A±39mA and at A5233A±84mA. 



intensity slope in the line is large and nearly linear at this point. Furthermore, 
this part of the line is well correlated with the two adjacent sampling points: 
the Pearson's r correlation coefficient has values of 0.50, 0.90, 0.90 and 0.50 for 
m = 9, 29, 102 and 177 and is 0.96 for b, m = 5250, 39. The coefficients and 
their uncertainties have been computed using the methods described by Isobe 
et al. (1990) and Feigelson and Babu (1992) wherein neither of the two quantities 
is treated as dependent or independent. The multiplicative factors to convert the 
magnetic fields to the A5233A ± 84mA scale are given in Table 1 and deviate 
from each other. This clearly indicates that A5233A does not provide a single 
reference field strength. The deviation of the scale factor for A5233A ± 29mA 
from unity is particularly significant in that these field measurements are well 
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Figure 4. Part a of this figure shows the spread of points on the scatter diagram (labeled 
Detronded) compared to the spread found from a standard magnetogram observed with the 
polarization modulating KDP analyzer inactive (labeled KDP off). Part b of this figure shows 
the ratio of the gaussian width of the detrended scatter diagrams for the five parts of the 
A5233A line shown as the symbols compared to a model that assumes that the variation in 
measured field is due to photon noise shown as the solid line. 



correlated with those from A5233A±84mA in spite of the fact that the hne wing 
profile is straight and neither is close enough to the line core for the flux tube 
field to shift the Zecman component past the sampling point. In the next section 
we develop an algebraic model for the saturation factor and then in section 5.2 
apply this model to observed profiles of the A5233A line to verify the above 
conclusion. 

One issue that arises later in section 5.1 is the cause of the spread in the 
scatter diagrams shown in Figures 1 to 3. Two possible causes are instrumental 
noise mostly due to photon statistics or an intrinsic deviation from a perfectly 
linear relationship between the different pairs of spectral sampling points. We 
distinguish between these two explanations by examining distribution functions 
(j){x) where 0(x) is the number of pixels having x between x and x + Ax and 
we use a normalization where 0(0) = 1. For x we use either x — B — Bm or 
x = Bkdp off - For _B ! s the value of Bf,f for each point is derived from 

A5250A±39mA ^ 

B ! ! using the trend line in Figure 3. For Bkdp off we use an obser- 

A5233 A.in84mr\. 

vation from the regular synoptic magnetogram program where for calibration 
and control purposes we periodically measure the magnetic field while the KDP 
modulator is not active. For the KDP off observation, the apparent magnetic 
field comes only from the system noise. Figure 4a shows this comparison. Clearly 
most of the scatter in Figure 3 comes from photon or system noise. 

A second part of the identification of the cause for the spread on the scatter 
diagrams is an estimation of the photon noise caused spread as a function of the 
spectral sampling separation AA. The deviation of V is the combined result of 
the deviations of the intensity measured for opposing states of the KDP analyzer 
with the form of the dependence being the same for all V^s. Each intensity has a 
fractional error that is inversely proportional to the number of photons collected. 
These intensity errors are converted into errors in V and subsequently into errors 
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in B. The factor converting V errors into B errors is inversely proportional to 
g{dl / dX) / 1 . The slopes and g values are known so the model for the relative 
errors only depends on the estimation of the relative number of photons measured 
for each spectral channel. All channels for A5233A use one set of fiber-optic image 
reformattors while the channel for A5250A uses an independent reformattor. The 
relative number of photons measured for all A5233A channels is just proportional 
to the relative line intensity at each spectral sampling point. Between these 
channels and the one for A5250A there is an uncertain relative transmission 
parameter. If we pick this single parameter to give the best fit then we can 
compare the spreads for all five spectral samples at A5233A to each other. This 
comparison is shown in Figure 4b. Clearly the model that the detrended magnetic 
field spreads is a consequence of photon noise is well established. 

3. The classic saturation effect 

The classical saturation effect comes from an interaction between an observing 
system and the fact that the observed pixels for the solar surface arc much 
larger than the magnetized portions of the solar atmosphere. Although the 
consequences of solar inhomogenieties have been intensively studied for many 
years (Solanki, 1993; Solanki, Inhester, and Schiissler, 2006), the shape of the 
A5233A line has led to the assumption that this line docs not suffer from the 
classical saturation effect. In view of the fact that we have found the observed 
magnetic fields derived from this line to depend on the spectral sampling, we 
feel it is important to quantitatively evaluate the classical saturation effect and 
definitively rule it out as a contributor to the spectral sampling dependence we 
find. 

The inhomogeneous character of the solar atmosphere means that the detected 
radiation is a combination with a fraction / from the flux tubes having a large 
Zeeman shift and a possibly different line shape and a fraction 1 — / coming from 
unmagnetized plasma with no Zeeman shift and a profile that may be closer to 
that from the quiet sun. The line shift amount AAh is related to the magnetic 
field strength by the standard formula: 

AAh = 4.67 X lQ-^^g^s\^B 

where A is in A and B is in gauss. (There is a notation conflict between this 
standard usage of AAj^and the quantity AA introduced in section 4.3 below. 
The usage with subscript H is restricted to this section and the later usages 
do not include this subscript.) Due to the non-linear shape of the line profile 
and the possible large shift in the flux tubes, it is not appropriate to calculate 
the observed change in the line intensity from f AXnidl /dX). In addition, the 
magnetogram observations do not measure the spectral line position for the 
two states of circular polarization since the movable stage does not translate 
at the rate at which the polarizations are switched. The stage servo drives the 
stage position towards the point at which the intensity averaged between the 
polarization states is equal for the spectral sampling ports on the blue and 
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red wings of the line. The position of the polarized line is then inferred from 
the Stokes V parameters at blue and red wing points on the line where the 
/ parameters are equal. Thus the operation of the system is then based on 
formulae originally discussed by Stenflo et al. (1984) with V being the difference 
in intensity for line profiles shifted by the Zeeman splitting amount. These can 
be represented using Taylor expansions for /(A ± AXh) up to third order. In 
the presence of an inhomogeneous atmosphere structure where the magnetized 
regions occupy a fraction / of the total area the formula for the V parameter 
can be written: 



Our interest is in the combination i?avo = /-Sft- The field strength and corre- 
sponding Zeeman shift AAj/p^ are taken to be consistent with a model external 
to our analysis and the deduced value of i?ave is then a function only of /. The 
third order relationship of interest for Save can be written 



where he quantity in front of /3(A) is just the standard formula for the magnetic 
field according to the weak field approximation. We denote the field derived from 
this term as -Bstokcs ^md use this uncorrected quantity as one of the primary ways 
of estimating the field strength. The quantity /3(A) is a saturation correction 
factor that depends on the line shape near the sampled point and the strength 
of the flux tube magnetic field which is given by 

Due to the difficulty in calculation of high order derivatives from observed data, 
we have replaced the first derivative of / by the Stokes V parameter. A similar 
replacement of the line slope by a function of V was recommended by Solanki 
(1993) for the study of flux tube field strength using the magnetic line ratio 
method. In our formula /3(A) depends only on the choice of flux tube magnetic 
field strength and not on the filling factor. The final formula based on V also has 
the advantage that it comes from the magnetized gas whereas the formula based 
on / includes radiation from both magnetized and unmagnetized gas. Due to the 
difficulty of taking the second derivatives, we only use this formula to estimate 
the size of the classical saturation factor and show that this effect cannot be 
responsible for the differences in magnetic field found from different parts of the 
line profile (see section 5.2.3 below). 

The above equations indicate the relationship between the circular polar- 
ization of the radiation and the magnetic field strength. However, the intensity 
emerging from individual flux tube which should be used in these equations is not 
available for observational determination under conditions of spatial resolution 
that prevail at the 150-foot tower telescope. The calibration observations carried 
out routinely prior to the magnetogram scans consist of spectral line scans in a 




V 



/3(A) 



avc — 



4.67 X 10-i3p3ffA2(a//9A) 
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quiet part of the solar atmosphere to yield /quiet (A). In principle the saturation 
factor /3(A) should be calculated from the intensity /ft (A) that emerges from 
a spatially unresolved flux tube. For reasons discussed by Frutiger and Solanki 
(2001) and as apparent from line profiles observed by our system and shown 
in Figure 6 below, /ft (A) ^ /quiet (A). Despite the line profile sensitivity to 
magnetic field, the tight correlations shown in Figures 1 and 2 demonstrate that 
the line shape effects are consistent over a large fraction of the solar surface so 
that the correction and filling factor effects in one area of the surface can be 
applied to other similar areas. 

4. Relationship between line profiles and magnetograms 
4.1. Context 

At each point of observation we wish to have an algorithm that returns a single 
value for the magnetic field strength even though we find a dependency of the 
deduced field strength on the profile sampling of the A5233A line. We verify below 
the expectation that the classical saturation effect does not explain this depen- 
dence so we must use other considerations to make this recommendation. As 
is well established (Solanki, 1993; Lites, 2000; Solanki, Inhester, and Schlissler, 
2006) the atmosphere consists of a weakly magnetized gas having one thermal 
stratification threaded by a varying density of highly magnetized flux tubes each 
of which has a different thermal stratification, the spectral line samples all of 
these structures with relative contribution factors that depend on position within 
the line. 

The strength of the magnetic field within the flux tube is generally thought to 
vary with altitude (Steiner, 1994) and the magnetic flux per pixel (the product 
of the field strength, the filling factor and the pixel area) may or may not be 
constant through the atmosphere since it is possible for field lines to close within 
the atmospheric zone. The radiation traversing the atmosphere containing these 
structures can be influenced by vertical and horizontal gradients that differ from 
the mean atmosphere. For example differently shifted Zecman components of the 
line might produce an effect similar to microturbulcncc for an emerging ray path 
over which the magnetic field strength is decreasing. These considerations lead 
to the fact that the shifted line profiles for the two states of circular polarization 
are generally not the same as the line profile in an unmagnetized region shifted 
by a constant amount and ftirther that they are not shifted by the same amount 
in all their parts. 

The line profile from an unmagnetized region gives us the residual intensity 
as a function of wavelength offset from line centre. An intensity difference such 
as that between the two states of circular polarization is then used to infer a 
line shift from the unmagnetized profile from an inverse where wavelength is 
considcrd to be a function of residual intensity. Although the shift of the flux 
tube line is large, the dilution with non-magnetized radiation makes the intensity 
difference also small so that the line shift can only be interpreted using the local 
non-magnetized line slope instead of the fully non-linear line shape information. 
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Table 2. The heights of formation, flux-tube field strengths and Zecman shift amounts. 
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122 
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92 
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22.5 


0.7 
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0.4 
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15.7 
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The difference between this interpretation and one that would be possible if we 
had adequate spatial resolution to isolate each flux tube is what we refer to as 
the classical saturation effect. 

4.2. Heights of formation 

Our task is to interpret the line profiles in terms of the flux-tube model and 
relate shifts in the line to a field strength appropriate to each part of the line 
and appropriate for each position on the solar disk. As a first step we have 
carried out a height of formation calculation for both A5233A and A5250A using 
the methods described by Caccin et al. (1977). We then deduced the standard 
optical depth from the atmosphere model of Vernazza, Avrett. and Loeser (1981) 
and inferred the magnetic field strength from the models of Frutiger and Solanki 
(2001) at the corresponding optical depth. The field strengths are given for our 
working points and center-to- limb positions in Table 2. These tables also give as 
AA// the line shift amount appropriate for the Zeeman splitting due to the field 
strengths at the working points. 

4.3. A line bisector definition of the magnetic field 

Solar surface magnetograms provided by the system at Mt. Wilson's 150- foot 
tower telescope are deduced from the longitudinal Zeeman effect starting with 
the measured intensity difference between two or more samples of a spectral line 
in the two states of circular polarization. As described by Howard et al. (1983) 
a double difference among the four measured intensities yields an indication of 
the opposing line shifts of the two states of circular polarization. The Babcock 
magnetograph includes a mechanism to track the Doppler shift of the observed 
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line and attempts to maintain equal average intensity for the spectral samples 
on the blue and red wings of the line. The mid-point between the blue and red 
equal intensity wavelengths is commonly referred to as a spectral line bisector 
point. The bisector points depend on position within the spectral line which can 
be defined either in terms of the ratio of the intensity to the continuum intensity 
(often called the residual intensity) or in the wavelength separation between the 
sampled blue and red wings of the line. Because the Babcock magnetograph 
defines a fixed wavelength separation, it is natural for us to adopt the bisector 
half-width AA as the parameter to define the position within the spectral line. 
(Recall that here AA without subscript H is not the Zeeman splitting used 
in section 3). This choice has the additional advantage that it is relatively 
unaffected by differences in the central residual intensity that we find to be 
common in comparing magnetized and unmagnetized regions. 

4.4. Magnetic fields from spectropolarimetry modelling 

Another way to interpret the polarimetry of spectral lines starts with the forward 
theory such as has been formulated by Unno (1956) and Landi degl'Innoccnti 
and Landi degl'Innocenti (1973). The polarized monochromatic intensity is then 
found from a solution to the transfer equation including the fact that the atoms 
in the magnetic field respond differently to the different states of light polar- 
ization. The resulting system of equations depends on the detailed stratification 
of the magnetic field and the physical state of the plasma - a highly complex 
combination that can only be managed with the help of approximations and 
simplifications. The fact that even the highest spatial resolution observations 
include contributions from magnetized and non-magnetized solar plasma means 
that a potentially very large number of free parameters is present in the task. It 
is especially unfortunate that spectra from individual flux tubes are not available 
so that the line shape from an individual flux tube must be recovered from the 
model rather than being observed directly. 

One form of solution is based on the use a Milne-Eddington approximation 
(Skumanich and Lites, 1987) and does not treat gradients in magnetic field 
strength. By including two or more spectral lines and invoking flux tube models 
with structural detail, Frutiger and Solanki (2001) have been able to estimate 
the thermal, dynamic and magnetic stratification of the flux tubes. Typically 
studies following this approach take the Stokes V polarization parameter as a 
function of wavelength near the spectral line as the primary constraining data. 
Works by Sigwarth et al. (1999) and Khomenko et al. (2005) show that a range 
of Stokes V proflles is encountered for high spatial and temporal resolution 
observations. As a consequence of this variability, it is likely that the underlying 
physical structure including the height dependence of the magnetic field is also 
variable. However, our interest is in the use of magnetograms for the purpose 
of establishing an inner boundary condition for global-scale magnetic models 
of the heliosphere. For this purpose, it is appropriate to use coarse spatial and 
temporal resolution which is less impacted by phenomena at or near the solar 
granulation scale. Indeed, the scatter diagrams of Figures 1 to 3 show that there 
is a tight correlation between the nominal magnetic fields derived from different 
parts of A5233A and between A5233A and A5250A. 
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4.5. Comparison to Advanced Stokes Polarimetry fields through MDI as an 
intermediary 

Although most studies using the techniques based on multi-line Stokes po- 
larimetry and inversions have concentrated on determination of the thermal 
and magnetic structure of the flux tubes, the work by Berger and Lites (2003) 
determined the relationship between the magnetic flux per unit area obtained 
with the Advanced Stokes Polarimeter (ASP) and the apparent field strength 
measured by the MDI magnetogram program. They found that the MDI fields 
are reduced to 64 ± 1.3% of the fields deduced from the ASP for the non-sunspot 
portions of an area centered around Active Region 8218 at S20 W22 on 13 May, 
1998 with a center-to-limb angle of 30°. For a similar position on the solar disk, 
the study by Tran et al. (2005) found that the MDI fields are reduced to 55% of 
the magnetic fields at A5250A measured by the MWO synoptic program after 
correction to the scale previously defined by A5233A. This comparison shows 
that the fields we previously deduced using the A5233A scale are 16% greater 
than those deduced using the multi-line spcctropolarimetry/invcrsion method. 

5. Line Profiles 
5.1. Observations 

To investigate the properties of the Stokes V profiles and deduced magnetic field 
strengths based on such profiles, 23 line scan sets were obtained over a period of 
several months during the declining phase of solar cycle 23. We use observations 
obtained between June and August 2007. At this stage of the solar cycle, the 
number of active regions available for study is limited. The observing system 
is that described by Ulrich et al. (1991) and utilized the 10-channel fiber-optic 
reformatter normally used for magnetograms in the A5896A line of sodium (see 
Ulrich et al. (2002) for the details of the spectral sampling). With this system 
we obtain spectral line profiles from a square selected by our entrance aperture 
which has sides of 12 arc-sec (other apertures are available but were not used 
in this study). This entrance aperture tracks solar rotation using the nominal 
differential rotation rate for the position observed. At the final spectrograph fo- 
cus the spectrum is sampled by one of four sets of fiber-optic image reformatters 
that are carried on two independently movable stages. A small portion of the 
solar spectrum including the target spectral line is scanned by moving its stage 
alternately redward and blueward over the solar spectrum repeatedly for a total 
duration of 20 minutes with each scan requiring 20 seconds. The second stage is 
left stationary and its output is used to remove sky transparency fluctuations. 
The length of each scan is lA and the spectral output from each fiber is sampled 
at an interval of 5.31mA (the spectral resolution is about 30mA). The temporal 
resolution was not used and all scans were summed to a single profile for each 
spectral pickup. Finally the ten profiles were shifted by the known flber-optic 
separation so that all are superposed for a single final profile. The offset between 
the fiber-optic inputs means that the portion of the spectrum scanned is also 
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Table 3. Line scan observations of A5233A 
Time is UT, Latitude, Central Meridian Distance (CMD) and 
Center to Limb Angle (CLA) are in degrees, the absolute value 
of the magnetic field \B\ is in gauss and has been determined ac- 
cording to the symmetrized Stokes V method below evaluated 
at i5A = 84mA. 



Set No. 


Date Time 


Lat. 


CMD 


1^1 


CLA 


022 


2007/07/13 


18:39 


-4.93 


5.56 


558 


10.7 


023 


2007/07/13 


18:52 


-2.89 


-0.35 


15.8 


7.2 


025 


2008/03/31 


19:01 


-11.32 


-54.49 


24.5 


54.0 


026 


2008/03/31 


19:24 


-2.43 


-53.17 


5.2 


53.1 


027 


2008/03/31 


19:56 


-10.83 


-70.45 


14.3 


69.60 


028 


2008/04/02 


18:08 


-5.26 


-80.60 


19.2 


80.02 


029 


2008/04/02 


18:30 


0.43 


-75.32 


2.3 


75.49 


030 




18:05 


-8.70 


-66.45 


68.6 


65.86 


031 


2008/04/05 


18:17 


-8.69 


-66.58 


63.2 


65.96 


039 


2008/06/14 


18:21 


-9.43 


2.48 


103 


10.82 


040 


2008/06/14 


18:51 


3.19 


-0.78 


16.2 


2.21 


041 


2008/06/15 


22:03 


-9.72 


-12.48 


8.4 


16.55 


042 


2008/06/15 


22:25 


-7.87 


-15.04 


132 


17.59 


043 


2008/06/15 


22:47 


-0.01 


-15.05 


12.6 


15.13 


044 


2008/06/21 


18:41 


-1.93 


-5.16 


161 


6.44 


045 


2008/06/21 


19:04 


-1.92 


3.34 


325 


5.18 


046 


2008/06/21 


19:25 


-1.93 


6.33 


3.8 


7.49 


047 


2008/06/22 


18:40 


-1.99 


-19.17 


314 


19.58 


048 


2008/06/22 


19:05 


-0.65 


-10.05 


174 


10.45 


049 


2008/06/22 


19:29 


1.22 


-0.60 


2.7 


1.02 


050 


2008/06/26 


19:06 


-1.20 


-62.56 


72.5 


62.64 


051 


2008/06/26 


19:28 


-1.89 


-71.48 


53.8 


71.56 


052 


2008/06/26 


19:50 


-1.33 


-57.32 


2.7 


57.41 



offset. The overall scan length was large enough that all ten inputs covered the 
range of interest for comparison to the magnetogram scatter diagram results. 
The circular polarization modulation was retained to permit the calculation of 
a single Stokes V profile. Following the convention described by Lites (2000) the 
blue shifted component is left circularly polarized in a positive magnetic field 
region. Although we retain knowledge of the polarity of each set of line profiles, 
we do not consider the sign of the field to be important in this application and 
present the Stokes V parameter and the derived magnetic fields as if all polarities 
are positive. Table 3 gives the parameters for these scans including the nominal 
magnetic field strength based on the line bisectors for a bisector width AA of 
±84mA. The line scan set numbers given in Table 3 are not sequential because 
a number of line scan sets for A5250A were also obtained during this time frame 
and these observations will be discussed in a future publication. 

The Stokes V profiles from the observations summarized in Table 3 are shown 
in Figure 5. The scale of the Stokes V plots is based on a continuum / of unity. 
The discussion in the following sections emphasizes the set of line scans from 
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Figure 5. This figure shows the Stokes V profiles for the 23 sets of fine profile scans. The scale 
of V is in percent relative to the continuum /. Each figure is labeled by the date of observation 
and sequence number allowing other parameters to be found in Table 3. The signal to noise 
of the polarimetry in these Stokes V profiles is approximately a factor of 10 greater than 
previously published cases and thus these measurements carry the analysis to field levels that 
are weaker than previously studied. In fact most of the variance in the V parameter as a 
function of wavelength and consequently also in the derived field strength comes from fringing 
effects introduced by our blocking filter. The estimated error in the derived field strength for 
the weaker cases is 0.1 gauss. The magnetic field as determined from the Stokes V parameter 
at it84mA is given on the figures as well as being included in Table 3. 



July 13, 2007 since that region was more highly magnetized than was the case 
for subsequent observations. The apparent noise for the low field strength Stokes 
V profiles comes from a combination of factors including photon noise and the 
formation of interference fringes by filters, windows and other plane-parallel 
optics in the system. We have mitigated these residual fringes by tilting the 
blocking filter but the effect is not eliminated by this technique. We estimate 
the uncertainty in the Stokes V profile over the spectral range needed for the 
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field determination is Vorror ~ 5.0 x 10^^. The line slope is 3 x lO^'^(mA)^^ at 
the points ±84mA from line center so that the Zeeman shift error due to this 
noise is about 1.7 x 10~^mA for an estimated uncertainty in the field strength 
of 0.1 gauss as observed by the channel at A5233A ± 84mA. We can also esti- 
mate the photon noise error in the deduced magnetic field using measurements 
from the normal synoptic program magnetograms. The integration time for each 
magnetogram pixel is 0.127 seconds while the integration time on each pixel in 
the spectral scans is very similar at 0.1 seconds. The final profile information 
is the result of combining scans from all 10 fiber-optic reformattor channels 
and from all 60 scans. Thus the relative photon noise for each line profile point 
compared to the synoptic program magnetograms should be reduced by a factor 
of (600/1.27)-'^/^ = 22. The line scanning process provides four times as many 
spectral points compared to the magnetogram leading to a further reduction of 
2 in the expected photon noise. According to the discussion of section 2, the 
gaussian width from the magnetograph distribution function is 4 gauss for the 
above spectral channel so that again we arrive at an estimate of 0.1 gauss for the 
spectral scans. For observations with the synoptic program channel at A5250A 
the equivalent error would be 0.025 gauss. These errors are substantially lower 
than has been previously reported and allow us to study the magnetic fields in 
regions that would normally be classified as unmagnetized. 

5.2. Interpretation 

For this section, we exclusively use the profile sequence number 022 of July 13, 
2007 to illustrate different ways of interpreting the profiles. Figure 6 compares 
the quiet and plage profiles. The two states of circular polarization are combined 
for the low magnetization profile but separated for the magnetized case. Because 
we are not concerned here with Doppler shifts, the low magnetization line has 
been shifted so its bisector position coincides with the average bisector position 
for the two magnetized components at a bisector half-width AA of 75 mA. This 
parameter measures the product of the filling factor and the fiux tube field 
strength and represents the magnetic flux emerging from the pixel. Not only do 
the two circularly polarized line profiles differ significantly from the quiet sun 
profile, they also differ from each other in such a way as to yield a magnetic flux 
which depends on position within the profile. The strength of this line and its 
high altitude of formation at line core make it difficult to attribute this effect 
to typical flux tube physics where thermal effects can produce line weakening 
or strengthening relative to the surrounding quiet atmosphere. Clearly then this 
spectral line docs not define a single magnetic fiux and our problem becomes 
that of deciding how to interpret the observations in such a way as to find the 
best estimate of the magnetic flux emerging from the solar atmosphere. 

5.2.1. Magnetic field estimated from the line bisector. 

The Stokes V profile is often used as the basis for this estimate. As an alternate 
way of thinking about the different parts of the spectral line it is worth con- 
sidering how a Babcock magnetograph might work. If we could simultaneously 
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track the spectral line in both states of circular polarization, the Babcock mag- 
netograph would evaluate the magnetic field on the basis of central wavelengths 
where the red and blue wing intensities are equal in each state of polarization. 
The wavelength difference between these two balance points would then define 
a magnetic field from the spectrally resolved line profiles. In general, a line 
bisector can be considered as either a function of the intensity relative to the 
continuum intensity or as a function of the wavelength separation between the 
equal intensity points. In view of the dependence of the line core intensity on 
the strength of the magnetic field, we find it advantageous to use bisector half- 
width AA as the defining parameter. Corresponding bisectors for the two states 
of circular polarization at a fixed but arbitrary value of AA are shown in Figure 
6. The separation between the two central wavelengths at this AA is indicated 
as (^Aj-jj from which a magnetic flux per unit area i3bi(AA) can be deduced. For 
the observation shown in Figure 6 we show the derived bisector field strength as 
a function of AA in figure 7. 

5.2.2. Magnetic field estimated from the Stokes V parameter. 

Figure 8 which gives both the Stokes V parameter and the slope of the quiet 
sun line profile. The scale factor for the Stokes V has been adjusted so that 
the curves agree near the line center. The magnetic field can be calculated from 
this scaling factor if we use the linear approximation. This gives a field strength 
based on the point-wise slope of the spectral line and the corresponding Stokes 
V value so that a magnetic field can be deduced from both wings of the spectral 
line. Wc refer to a magnetic field deduced this way as _Bstokos('5A). These fields 
along with the values of i?Bi(AA) are shown in Figure 9. Note for this plot 
we use 5X = ±AA. While -Bstokos ((^A) can be calculated for both red and blue 
wings of the line, i?Bi(AA) applies to both wings equally and has been plotted 
symmetrically. Although the Stokes V parameter gives independent magnetic 
fields in the blue and red wings of the line, a Babcock magnetograph responds 
to the average magnetic field from the two line wings. Consequently, apart from 
figures 8 and 9, the magnetic field derived from the Stokes V parameter has been 
made symmetric by averaging the blue and red wing field values. The comparison 
between i?stokes (<5A) and i3Bi(AA) indicates that the pointwise slope leads to a 
smaller deduced field strength than the bisector by 25% near the line core, a 
larger field in the mid part of the line wing then good agreement in the outer 
line wings. We do not have a definitive interpretation of this result but point out 
that it could be a consequence of microturbulcnt broadening of the line core in 
the flux tube. The influence of variable microturbulcnt line broadening on the 
interpretation of spectropolarimetry of A5250A has been discussed recently by 
Socas- Navarro et al. (2008). 

5.2.3. Classical saturation factor estimation. 

Because the line wings are nearly linear over much of the line, the classical 
saturation effect does not influence the measurement. The classical saturation 
factor derived from the Stokes V parameter is shown in Figure 10. In order to 
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Figure 6. This figure compares the line profiles from a quiet sun region to the two circularly 
polarized components from an active region which was magnetized with a negative polarity 
field. The quiet sun profile is shown as the solid line while the LCP and RCP polarized profiles 
are shown as the short dashed and long dashed lines respectively. For this figure the polarized 
profiles are normalized to the continuum being unity instead of 0.5 as is required for the 
calculation of the Stokes parameters. Also shown with the same dashed/undashed indication 
are the line bisectors for each of the three profiles. The parameter AA defining the position 
on the profile through the half-width value is shown for the two circularly polarized profiles. 
Finally this figure also shows the parameter (5Abi which gives the Zeeman splitting of the profile 
at the selected value of AA. The magnetic field at this position on the profile is found from 
<5Abi using the standard formula for the Zeeman effect and for this active region according to 
Table 3 the field strength was 558 gauss. 

illustrate /3 for a case of interest, a flux tube field of 1500 gauss was used for 
this calculation although the appropriate field strength for the line core is closer 
to 400 gauss according to Table 2. Even with this higher than appropriate field 
assumption, the value of [3 deviates from unity by only 10%. The A5233A has 
nine Zeeman components. As a simple numerical test, the unmagnetized profile 
was shifted and weighted by these nine amounts and a simulated profile was 
then produced. This numerical test gave essentially the same result as we find 
by shifting the line an amount equal to the weighted mean of the Zeeman factors. 
We conclude from this simulation that the saturation factor is not substantially 
influenced by the presence of the nine sub-components. These parts of the line 
may however play a role in the line transfer process in a manner that is not 
captured by our simple approach; the multiple Zeeman components may produce 
the equivalent of an excess of microturbulence. We conclude that the dependence 



Magnetic_Scale.tex; 12/12/2008; 5:39; p. 19 



Ulrich et al. 



500 ^ 

400 

5 300 
CO 

O 

CD 200 



1 00 U 

r ' 

50 100 150 200 250 300 

AX (mA) 

Figure 7. The magnetic field calculated from the difference in wavelength of the line bisectors 
of the two states of circular polarization. The position within the line of the bisector is defined 
by the bisector half width AA. Note that this field differs from that given in Table 3 wherein 
the fields from the Stokes V parameter are given. 
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Figure 8. A comparison of the slope of the spectral line in the non-magnetized comparison 
region to the Stokes V parameter from the nearby magnetized region. The Stokes V parameter 
has been scaled to have the same shape near the line center. In the positively magnetized region 
the observed Stokes V parameter reached a maximum of 3.7% at -72 mA and a minimum of 
-3.1% at 4-70 mA. 

of the deduced field strength on position within the spectral line for A5233A is 
not consequence of the classical saturation effect even for points near the line 
core. 

The P parameter can also be derived for the A5250A line. In this case the 
greater curvature of the line profile leads to /3 = 1.6 at our working point of 
±39mA. The factor is only 35% of the slope of the relationship indicated by 
figure 3 confirming that thermodynamic effects play a dominant role in the 
modification of the apparent magnetic field strength. 
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Figure 9. This figure shows the ratio of the Stokes V parameter divided by the spectral line 
slope. The left scale gives the result in units of mA while the right scale gives the inferred 
magnetic field strength in gauss. In this figure, the asymmetric V profile is used. Also shown 
is the magnetic field deduced from the line bisector separation using the bisector width AA 
as being equal to the distance from line center SX. This function is plotted for both +AA and 
— AA since it is symmetric as derived. 
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Figure 10. This figure shows the classical saturation factor /3 for the line A5233. The flux 
tube field was taken to be 1500 gauss as is consistent with points in the wings of the line. Near 
line core where the field strength is about 400 gauss, the value of f3 is much closer to unity. 



6. Magnetic fields referenced to A5233A 



At each point on the solar surface we now have a range of estimates for the 
magnetic field: T^^f^^'^^Bb.m with 6, m = 5233, 9, 5233, 29, 5233, 102, 5233, 177 
and 5250,39, i3Bi(AA) and i?stokcs(<5A). We should choose the estimate which 
gives the product of the field strength and filling factor appropriate to the pixel 
area contribution to the magnetic flux above the solar atmosphere. For this 
purpose, the portion of the line nearest the line center would appear to be the 
best choice since this part of the line is highest in the atmosphere. However, as 
we show in this section, the choice has to be tempered by the fact that the line 
core shows an excess range of variability and is vulnerable to line broadening 
uncertainties from microturbulence. 

For portions of the solar surface outside of sunspots the apparent magnetic 
field strength is not greater than about 500 gauss. For Zceman shifts from fields 
of this strength the linear approximation is applicable and the field deduced from 
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Figure 11. This figure shows comparisons between the field strength dependence on position 
in the spectral fines and the slopes of the correlations given in Table 1 . Each panel corresponds 
to a line scan set from Table 3 as indicated on the panel. The solid lines are for -Bstokes(^-^) 
and the dashed lines are for i?Bi(AA). The Stokes V profile from Figure 9 has been made 
symmetric by averaging the blue and red wing results since we are applying these results to 
magnetograph observations where the blue and red wings are averaged by the instrument. The 
field strengths in gauss for both the bisector separation and the Stokes V ratio divided by the 
slope are given. The relative correlation slopes have been scaled to correspond to the field from 
the Stokes V ratio. The x symbols are from Table 2 while the + symbol is taken from the 
report by Uhich (1992). 



a magnetograph observation depends on the spectral hne slope at the point of 
the line where the spectrograph sampling is selected. This means that the magne- 
togram field is consistent with what we derived above from the ratio of the Stokes 
V parameter to the local spectral line slope - i.e. the magnetogram magnetic 
field should be consistent with the Stokes V magnetic field interpretation of 
the spectral line rather than the bisector interpretation of the spectral line even 
though the the Babcock magnetograph follows the average spectral line bisector. 
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Thus, it is reasonable to expect that the observed magnetogram correlations 
shown in figures 1 to 3 and Table 1 between different deduced fields should 
be consistent with the Stokes V deduced fields. The line profile observations 
come from a relatively small number of selected regions on the solar surface 
that may or may not be representative of a general case. However, since many 
of the slopes from different spectral sampling configurations are well correlated, 
these particular observed profiles must be representative of the general case since 
otherwise, there would be a wider scatter in the correlation diagrams. Thus the 
implied field values from different parts of the line profile can be normalized so 
that the Stokes V field strengths and bisector field strengths give scale factors 
that can be used in bringing the different deductions to a common value. 

The 23 sets of line scans have been reduced to yield the three estimates of 
the field strength indicated above. The value of i?5233,84 was chosen to pro- 
vide a set of field strengths roughly similar to the other two estimates. The 
value of ?75 2 33'84 ^''^^ interpolated in sin(p) to be appropriate to the position of 
each line scan field strength function. Selected line scan results illustrating this 
comparison are shown in Figure 11. The X 's indicate the rescaled correlation 
coefficients obtained from the magnctograms, the solid and dashed lines give the 
symmetrized fields from Figure 9. In addition to the cases of 775233 §4 from Table 
1 we have included the implied slope at A5233A±45mA by combining the slope 
in Figure 3 with the correlations from the observations of 1991. This point is 
shown as the + symbol. 

The choice of cases to present in Figure 11 was made to illustrate four points: 
1) the weaker field cases are erratic, 2) the bisector fields show less dependence 
on position in the line than is the case for the Stokes V field, 3) for the stronger 
field cases the slope dependence on spectral sampling is similar to the Stokes 
V dependence on 6X and 4) there docs not seem to be a dependence of the 
results on the ccnter-to-limb angle. For field strengths above about 40 gauss, the 
patterns for the cases shown in the top four plots are typical. For the weaker 
fields there is a considerable range of behavior with that for the bottom left case 
being an extreme. The implied line core field for this case is only 20% of that 
implied at AA = 84mA. Other cases of weak field actually have the core field 
greater than that at AA — 84mA. In contrast the higher field cases show the 
pattern of the bisector field being less dependent on AA than the Stokes V field 
is on SX for portions of the line nearer the line core. 

Since the Stokes V field Bstokcs gives a good reproduction of the correlation 
coefficients from the scatter diagrams and since the core parts of the bisector 
field i?Bi(AA) are less sensitive to AA than -Bstokcs('5A) is to SX, we adopt a 
model that the most appropriate field to use in determining the magnetic flux 
above the solar atmosphere is i?Bi(AA) where AA approaches zero. We believe 
that correlation coefficient field as well as the Stokes V field near the line core 
are influenced by extra line broadening in the magnetized flux tube. Since we 
are unable to observe this line profile we must use the profile from a standard 
quiet region which has a steeper slope near the line core so that the apparent 
falloff in field strength is due to division by too large a number. It is possible 
to estimate some properties of the spectral line profile in the magnetized flux 
tube by integrating the Stokes V parameter. This gives a shape near the line 
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Figure 12. This figure shows the factor y which converts the slope ^^250 39 correction 
factor that will yield our recommended magnetic field. The x symbols represent line profiles 
observed with center-to-limb angles less the 45° while the + symbols are for profiles with 
center-to-limb angles greater than 45°. One point with B = 558G and y = 0.689 is not shown 
to avoid compressing the abscissa scale. This point has a center-to-limb angle of 10° . 



core which is indeed broader compared to that of the unmagnetized hne in the 
case of the observations of July 13, 2007. The bisector field does not depend 
on the line slope but only on the requirement of equal intensity on the red and 
blue wings so that its falloff toward the line core must be a consequence of some 
other effect such as flux cancelation. To combine the results of the magnctogram 
scatter diagrams and the line profile data using our adopted model, we identify 
the line profile quantity Sgtokcs (84mA) as corresponding to -65233,84 used in 
the scatter diagram and we identify the line profile quantity i3Bi(29mA) as the 
recommended field. The ratio: 

BBi(29mA) 



Sstokcs(84mA) 

is then a factor which will convert t75 25o'39 into the final correction factor for our 
observed magnetograms using A5250A ± 39mA in the MWO synoptic program. 
We show the derived values of y in Figure 12 as a function of i3stokos(84mA). 

It is noteworthy that for areas with fields i3stokcs(84mA) larger than about 30 
gauss the values of y are quite stable and do not depend on the center-to-limb 
angle. The variations in y for smaller fields are not due to error of measurement 
of the fields but rather are a consequence of some structural property of the 
solar surface. For our purpose, we do not need to understand the cause for this 
behavior and can simply adopt the value of y to be 0.73 ± 0.015 which is the 
average and error of the mean of all the cases with i?stokes(84mA) > 30 gauss. 
Examination of Figure 5 reveals that the raw V profiles for the weaker field 
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Figure 13. This figure shows the final scale factor <5~^ which converts the observed mag- 
netic field for A5250A±39mA into the recommended magnetic field strength that gives our 
best estimate of the magnetic fiux from each pixel emerging into the region above the solar 
atmosphere. The x symbols with the error bars give the values of f)5233'84 with their formal 
uncertainties multiplied by the average value of y of 0.73. The solid line is the fit to sin^(p) 
with the fitting coefficients given above the line. In a similar fashion the filled triangles give 
the points from Ulrich (1992) with dashed line showing the fit to these points and the fitting 
coefficients in use by Wang and Sheeley (1995). 

cases are much less regular than the stronger field cases. The flux tube structure 
implies that the weaker field areas have a small number of flux tubes. The V 
irregularities may be a consequence of statistical fluctuations around a fairly 
stable mean structure. Theoretical studies of the A5233A line would be helpful 
in understanding the properties of this line in the magnetized and unmagnetized 
solar atmosphere as well as providing a separate indication of the relationship 
between the various observed magnetic held strength indicators and the magnetic 
flux emerging into the heliosphere. 

Based on the adoption of the i3Bi(29mA) as the recommended field value, we 
can use the magnetogram correlation diagrams to obtain a new scale factor for 
the Mt. Wilson 150-foot tower telescope synoptic program which utilizes the Fel 
line at A5250A. For notation, we retain the traditional quantity first defined 
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by Howard and Stenflo (1972) wherein the recommended field Bjiec is: 

^Roc = <^ ^^5250, 39 ■ 

The previous correlation analysis by Ulrich (1992) as represented by Wang and 
Sheeley (1995) was: 

5-'^ =4.5-2.5sin2(/9) . 
The present analysis is shown in Figure 13 and gives: 

S'' = ^5250,39 = 4.15 - 2.82sin2(p) . 

Note that this correction factor is simply to take into acount line formation 
processes and does not include effects related to the geometry of the field. 
In particular an additional sec{p) correction is required if one adopts a model 
wherein the magnetic field is radial. 

7. Discussion 

Our principal result is given above and provides a small correction to the proce- 
dure previously applied to the Mt. Wilson Observatory magnetogram observa- 
tions. The previous description of this correction factor as a saturation correction 
is false. The various interpretations on the polarimctry of the spectral lines 
used in our observations differ from one another due to processes involved with 
transfer of radiation through an inhomogencous, magnetized atmosphere. A full 
analysis of this problem is beyond the scope of the present investigation but 
is obviously needed to fully understand the changes in shape of the polarized 
spectral lines. We base our recommended correction factor on the fact that the 
magnetic field deduced from the spectral line bisector is less variable than the 
field deduced from the monochromatic Stokes V parameter. We further adopt 
the magnetic field implied by the properties of the spectral line near its core 
on the grounds that this part of the line is closest to the outer edge of the 
solar atmosphere. We find that the line profile properties are most stable in 
regions where the field strength is significantly larger than that found in the 
completely quiet sun. The line profiles for the weakly magnetized regions are well 
determined by our system but they show substantial deviations from symmetry. 
We suspect that this is a consequence of the reduced effect of adjacent flux tubes 
as a constraint on the flux tube structure in cases where the field is so weak that 
there are no adjacent flux tubes. The unconstrained flux tubes can then adopt 
a wider range of structure including perhaps cases where the field lines return 
locally to the solar interior rather than extending out of the solar atmosphere. 
The fact that the error of measurement of the fields is well below 1 gauss for our 
line profile system means that the measured profiles can be used to constrain 
models of fiux tubes in regions where the field strength is very low. 

Our determination that the differing deduced magnetic field strength for 
differing spectral sampling configurations has a general potential impact on 
the measurement of magnetic fields with other observing systems. All programs 
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provide maps which give a single magnetic field strength for each observed point. 
Had the adjustment of the Mt. Wilson fields been a consequence of a saturation 
effect unique to this system, the potential implications would also have been 
limited to the Mt. Wilson system. However, our determination that the effect 
comes from a spectral transfer effect means that other observing programs need 
to determine how to relate their field measurements to measurements made 
by other systems. Even within a single observing programs, multiple deduced 
magnetic fields could occur depending on the exact specification of the algorithm 
bringing the spectropolarimetric data to a deduced magnetic field strength. Ex- 
amples where a direct comparison has been made on a pixel by pixel basis include 
Tran et al. (2005) and Demidov et al. (2008). We are also in a position to update 
discussion from 4.5 which compares our scale factor to the deductions of Berger 
and Lites (2003). Those authors concluded from Advanced Stokes Polarimeter 
data that the MDI magnetic fields are 0.64 ± 0.013 times the correct value. 
After combining the factor from Tran et al. (2005) to the new corrected A5233A 
scale we find that the MDI magnetic fields are 0.619 ±0.018 times the recom- 
mended value at the center-to- limb distance used by Berger and Lites (2003). 
This agreement within errors of measurement gives support to our approach in 
the calibration of the A5250A magnetograms. 
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